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Abstract 

We present some arguments showing spectrum doubling of matrix models in the 
limit ./V — > oo which is connected with fermionic determinant behaviour. The 
problems are similar to ones encountered in the lattice gauge theories with chi- 
ral fermions. One may discuss the "physical meaning" of the doubling states or 
ways to eliminate them. We briefly consider both situations. 
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1 Introduction 



A tool for nonperturbative study of superstring theory (M-theory) proved to 
be matrix models, in special BFSS (M(atrix) theory) and IKKT (IIB matrix 
theory) ones, [1,2]. These models are formulated in terms of N x N Hermitian 
matrices, where N — > oo. 

For a finite (but large) N the last model, [1] plays the role of regularised 
"second quantised" IIB superstring in Schild formulation [3] , while the former 
one [2] is a regular isat ion of D = 11 quantum super- membrane. 
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Even for finite N these models have problems with definiteness due to diver- 
gence of the partition function arising from integration along flat directions 
or vacua. Due to this also the limit N — » oo is bad defined. However, if not to 
try to consider the matrix model as a whole, but to consider a "perturbative 
sector" connected with fluctuations around a particular vacuum solution one 
may have a well defined system. 

The models obtained this way depend on the chosen vacuum solution, in 
particular, various compactifications of IIB matrix model on noncommutative 
tori were shown to yield super Yang-Mills (sYM) models on these tori, as well 
as BFSS model, [4]. 

One can show that in fact Yang-Mills type models describe all regular fluc- 
tuations around a vacuum solution of IKKT matrix model, [5]. Under regular 
fluctuations we understand those which have bounded value of momentum 
and position operators. In practice e.g. numerical computations this restriction 
can be implemented by additional regularisation which makes the undesirable 
modes decouple. 

Once there is a close correspondence between matrix models and the sYM 
models there appears the following problem. Ten dimensional sYM model is 
known to be anomalous and the anomaly seems to exist also in the noncommu- 
tative case, [6,7]. From the other hand there is no visible source for anomaly 
in the matrix model. For any finite N, the matrix model is finite manifestly 
gauge and Lorentz invariant as well as supersymmetric. Naively, these proper- 
ties must hold also in the limit iV — > oo, contradicting the anomaly of D — 10 
sYM. 

The explanation of this discrepancy may reside just in the contribution of the 
singular configurations. In the actual paper we show that the matrix model 
fails to reproduce the anomaly in the limit iV — > oo due to the spectrum 
doubling in the fermionic sector of the matrix model. 

In fact the actual situation is not new. A class of models is known to suffer 
from the spectrum doubling, [8,9]. In particular, for lattice gauge theory, it 
is known that the the naive discretisation of chiral fermionic action leads to 
doubling of the fermionic spectrum. The last results in appearance for each 
mode another mode(s) carrying the opposite chirality which compensate the 
parity odd contribution of fermions, [10]. In fact, there is a no-go theorem due 
to Nielsen and Ninomiya [11,12], which states that the doubling cannot be 
avoided unless the gauge, Lorentz or other relevant symmetry is destroyed in 
the continuum limit. While the gauge symmetry plays the crucial role in the 
consistency of the model, breaking of the Lorentz symmetry is not danger- 
ous for the model and maybe even desirable as a mechanism of spontaneous 
breaking of the ten dimensional Lorentz group to a lower dimensional one, [13]. 
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The plan of the paper is as follows. First we introduce the description of the 
fermionic sector of IIB matrix model for finite N, after that we explicitly find 
the doubling states for free fermionic fluctuations (quadratic approximation), 
and discuss the issue for the interacting case as well as possibility to eliminate 
the doubling. 



2 Finite N Matrix Model 

The IIB matrix model is described by the action, 

S = ~tr QK, A,] 2 + ^L4 M , V#) , (1) 

where and ip are ten dimensional vector and Majorana-Weyl spinor Her- 
mitian N x N matrices. 

An important class of configurations which tend to satisfy equations of motion 
in the limit N — > oo, are BPS ones given by = p M , where Hermitian 
matrices p M = satisfy, 

\Pn,Pu]=iB IJU/ , (2) 

here B^ u is proportional to the unity matrix B^ u = ■ I. 

Although, such a set of matrices does not exist for finite N, it can be approx- 
imated by a sequence of matrices converging to (2) in the sense of operator 
norm on the Hilbert space of smooth finite functions (vectors), [5]. 

Before searching such an approximation, consider a Lorentz transformation 
which brings matrix B^ u to the canonical form, having 2x2 antisymmetric 
diagonal blocks with ±^ entries. In this case the set of matrices p^ is split in 
pairs (pi, (f ), i = 1, . . . , D/2, where Pi and q l are canonically conjugate, 

\Pi,Pj] = [q\q j ] =o, (3) 

\p t ,qi} = -m t 5l (4) 

Eqs. (3,4) give the D/2 dimensional Heisenberg algebra. It is known that 
the Heisenberg algebra can be represented e.g. in terms of square integrable 
functions defined on on the spectrum of operators q l . 

Finite iV analog of the Heisenberg algebra (4) is given by the position and the 
Hermitian (symmetric) lattice derivative operator on a compact rectangular 
periodic lattice Y. Let |n) be eigenvector of q % with eigenvalues Lj sin(27r/iVj), 
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g*|n) = L i sin(27rn i /^i)|n), 



(5) 



then operators pi act on this basis as follows, 



Pi\n) 



itii 

2a,i 



(\n + e t ) 



n - d)) , 



(6) 



where is the unity lattice vector along i-th link, and Oj = -jj^- The choice 
of q which we use differs from one usually used in lattice model by redefinition 
of q — > Lsm2nq/L. It is as good as the former one, but beyond this it treats 
p and q in a symmetric manner. It is not difficult to see that in the basis of p 
eigenvectors eqs. (5) and (6) keep the same form with p and q interchanged. 
It also introduces explicitly the periodic boundary conditions, i.e. the identifi- 
cation \n + iVjej) ~ \n). Due to the symmetry between p and q the continuum 
limit is achieved after the "UV cutoff" removing: a« — > 0, combined with the 
"IR cutoff" removing or "decompactification" : L; — > oo. 

The commutator of and g J looks as follows, 



where we have introduced the notation I^)(jV) for the following matrix, 



As it is not difficult to see the equations of motion are not satisfied by such 
background. It still can be shown that the N = oo solution (2) can not be ap- 
proximated by solutions to equations of motion at a finite value of N, because 
at finite iV the only solutions are those with zero commutator, [A^ A v \ = 0. 

For (sequences of) vectors |/) = J2n fn\n), on which p and q remain bounded 
as N goes to infinity, 



x 




i^I a) (A0|n>, (7) 




(f\(p 2 + q 2 )\f)<C, 



(9) 
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where C does not depend neither on N nor on L (or a), operator %)(-/V) 
approaches the unity one, 

%>(W>« 

2^ 2 U n+ei COS AT + /n ~ e * COS AT ) l n ) ^ I/)' ( 10 ) 

since in this case n? <C iVj, and f n ± ei — f n — O^N^ 1 ). 

The last equation means that the operators preserving the property (9) tend 
to commute with I^(N) as iV approaches the infinity. 



3 Fermionic Contribution 



We are ready now to proceed to the analysis of the fermionic contribution to 
the partition function of the model with action (1). Integration over fermionic 
matrices results in the Pfaffian of the fermionic operator, 

Z(A) = J (tye-* (11) 

where, Sf stands for the fermionic part of the action (1). 

For finite N consider the bosonic background given by matrices pi, q l from eqs. 
(5,6). An arbitrary Hermitian matrix fluctuation around the given background 
is = p^ + ga^. The fermionic part Sf of the action in this case looks as 
follows, 

Si = -^tr^H(p M + ^),V], (12) 
where we rescaled if) — > gtp. 

The free (a M = 0) part of the fermionic action can be written in the represen- 
tation of \n). It looks as follows, 

S{ = [{-ih i /(2a))i) ntin r('4} m+eun - i) m - eun - i> m , n+ei + w„_ e J 

m,n,i 

r 7 S ( ■ 2nmi 2 ^ ni \ , \ 

+L i ip ntm T i I sin — sin I ^ m ,„ I , (13) 



where, 



ipn, m = {n\ij)\m), (14) 
i> m ,n = {m\^\n), (15) 
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and P, Ti are Dirac matrices, 



(16) 



Although, the action (13) differs from the naive lattice fermionic action, it 
shares many common features with it. 

As for naive lattice fermions, they are described by the action [10], 

Snaive = ^ £ (VVf-e; ~ VV-eJ, (17) 

^ n 

where for shortening notations we put lattice spacing a to unity. 

A well-known fact is that the actual model suffers from the fermionic spec- 
trum doubling. The last manifests in the fact that for each chiral fermionic 
state there is always another one present in the spectrum with the opposite 
chirality but with other quantum numbers coinciding with the original state. 
This phenomenon can be described by introducing some discrete symmetry 
which relates these states. The states obtained by action of this symmetry 
are called doublers. It is clear that if such a symmetry exists it completely 
destroys the chiral asymmetry.^ 

In the case of action (17), such a symmetry indeed exists and its generators 
in a even dimension D look as follows [14], 

T a = T {D+1) r a (-l) n °>, (18) 

where Tr D+1 \ is the D-dimensional analog of the Dirac 75-matrix (D is even), 
= cdT 1 !" 2 • • • T D , e is chosen to be either i or 1 in order to make F(d+i) 
Hermitian. 

Finding the order of the discrete group generated by T a , one finds that the 
number of doubling states is 2 D — 1. 

Now, let us return back to to the matrix model given by eq. (13) and try to 
find a similar symmetry in this case. 

One can check that the action (13) is left invariant by the following symmetry, 

^ m , n - ir n r,(-ir l - m > m , n . (19) 

or in the matrix form, 

V -»• Tt ■ V = iTnTiUr^Ui, (19') 

3 The absence of an explicit symmetry of such kind, however, does not prove nec- 
essarily, the absence of doubling. 
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where unitary matrix Ui is given by [/, = (— 1)^ arcsm (g7Li), 

Indeed, factor TnTj commute with all f j and T- 7 for j ^ i while factors 
(— l) nl_ml are the same for both ip and ip, and, therefore cancel. In the re- 
maining term TnTi anticommutes with I\, but the extra minus sign is com- 
pensated by the variation of the factor (_i) n '- m \ Thus, all terms in the action 
(13) remains invariant under the transformation (19). 

Interchange p <-> q gives the remaining symmetries, 

i) - % ■ i\> = iTuTiUr^Ui, (20) 
where Ui is the unitary transformation, Ui = (— 1) 27 arcsm (P ai /^). 

Summarising one has the action (1) invariant with respect to discrete symme- 
try generated by T M , 



where C/„ satisfy, 



T,-^ = iT n T,U^U„ (21) 



U^PuUn = (1 - 2^K- (22) 



As in the case with naive lattice fermions these transformations act in such 
a way that in the continuum limit the states become 2 D -fold degenerate with 
half of that for each chirality. In particular eq. (20) means that for each matrix 
state of given chirality which connects m and n there are 2 D / 2 states of different 
chiralities connecting ^e— m with ^e^— n, ^ej+^ej+m with ^-ei+^-ej+n, 
i 7^ j, and so on. Eq. (19) have the same interpretation in the "momentum 
space" spanned by the eigenvectors of p. 

It is clear now that if one wants to compute the gauge anomaly, one will have 
contributions from doublers of both chiralities which cancel each other. 

So far, we considered the interaction free part of the fermionic action. Presence 
of interaction at least in the framework of perturbation theory does not change 
the situation, in this case one has in the continuum limit an interacting 2 D - 
plet instead of a free one. Let us note, that this analysis may not remain true 
beyond the perturbation theory, as for strong field a M the interaction part 
dominates and the symmetry (21) of the free part does not play in this case 
such an important role. Unlike the usual lattice models the study of doubling in 
nonperturbative regime is too complicate, but the perturbative considerations 
are enough to doubt the result of naive continuum limit. 

One may ask, what happens in this case to the supersymmetry of the original 
matrix model given by the action (1)? 
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In fact, in spite of the fermionic doubling, the supersymmetry still exists as 
it is valid for any matrix configuration and is irrelevant to the chosen repre- 
sentation. The apparent paradox with doubled number of fermions is solved if 
one see that the number of bosonic "degrees of freedom" is also doubled. The 
bosonic doublers are gauge equivalent configurations and are related by gauge 
transformations a M — > U~ 1 a fl U u . 



4 Doubling States Removing 

The results of the previous section show that matrix model (1) fails to repro- 
duce a chiral continuum model in the limit iV — > oo. Thus, in particular, one 
can not obtain from it the noncommutative SYM model in this limit. 

Such a situation can be interpreted either as a presence of finite N artefact 
which does not decouple in the limit N — > oo, and must be removed by ad- 
ditional effort like in traditional lattice models, or as an indication that the 
model possesses nontrivial symmetries. In traditional lattice models the dou- 
bling states should be removed since the doubling contradict the continuum 
"phenomenology", in special the chirality properties of the model. 

In the case of matrix models which pretend to describe M-theory the situa- 
tion is different. As it is known the M-theory unites perturbative models with 
different field content. In particular IIA models contain states with both chi- 
ralities while in IIB models there are only ones with definite chirality. The 
duality symmetry which must relate these models should contain a mecha- 
nism which flips the chirality of states. As we see, such a mechanism exists on 
compact noncommutative spaces and is given by doubling. 

In spite of this perspective for doubling states to describe the physical reality 
there exists, however, possibility to remove them in order to get a chiral model 
in the continuum limit. Consider briefly the ways one can do this. As the 
problem is a "lattice" one, we can look for specific lattice solutions. In the 
lattice case the doubling is cured by addition of a Wilson term to the naive 
lattice action [10]. The problem is that there is no gauge invariant Wilson 
term which could be added to our "naive" fermionic action (12), as there is 
no gauge and Lorentz invariant fermionic mass term in the model. 

One can, however, write down Wilson terms preserving either of two symme- 
tries. 

A possible gauge non-invariant Wilsonian prescription is given by addition to 
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the naive action (12) of the following term, 

AS w , g augc = -^rfjF^ri] + [p^Mp^H] + \p { +\ W^W, (23) 

where r\ is U (iV)-singlet Majorana-Weyl spinor matrix, and p^ are respec- 
tively forward and backward one-step scaled lattice derivatives, 

pl ±) \n)=±d-(\n±e i )-\n)) (24) 

V «i 

q {±)i \n) = i/^e ±2 ™'/^ \n). (25) 

Due to the term (23) the states with large phase of q and p, (i.e. with n l ~ iVj 
and fcj ~ ^) acquire large masses and decouple in the limit N — > oo, as it 
happens in the case with usual Wilson term. 

Another possibility is given by that in contrast to low- dimensional field theory 
models where the Lorentz invariance is obligatory, in the Matrix model it is 
less important. Moreover, its breaking to lower dimensional symmetries is 
desirable if one wants to describe a four dimensional theory in low energy 
limit, as it was proposed in the Ref. [13]. It is not difficult to construct Wilson 
term which breaks, say Lorentz group SO(10) down to SO(9), but preserves 
the gauge symmetry. It looks as follows, 

AS w ,Lorc„t z = ~\ tr[pW, 4>}r%~\ # (26) 

where T 9 is the 9-th ten dimensional Dirac gamma matrix. 

Under this choice modes with large n and k also acquires large masses in the 
limit N — > oo, but it produces terms which are not invariant with respect to 
rotations involving the 9-th axis. The gauge symmetry here remains intact. 
The last should not appear strange, because giving up a part of Lorentz in- 
variance in an anomalous model allows one to cancel anomaly. Thus, in the 
simplest case of Abelian gauge anomaly in D = 4 one can cancel the anomaly 
e^^F^Fxrj, where //, u, ■ ■ ■ = 0,1,2,3 by addition a local counterterm 
~ A Q e l i k AidjA k , with = 1,2,3 which is not invariant with respect to 

Lorentz boosts. 



5 Discussions 

We have shown that perturbative fluctuations in IIB matrix model at finite N 
exhibit a phenomenon similar to one in lattice gauge theories with fermions, 
consisting in doubling of the fermionic spectrum. 
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We considered a simple example of a background given by lattice shift/mo- 
mentum shift operators. Basing on analogy with lattice model, we conjecture 
that this is a universal feature appearing for arbitrary choice of background 
configuration which is Hermitian, nondegenerate, etc., and can not be elim- 
inated without breaking either Lorentz or gauge invariance. The actual results 
are straightforwardly translated to the case of the genuine finite N vacua con- 
sidered in [5]. 

We give prescriptions for the elimination of the doubling states in the limit 
N —>■ oo, but preserving either Lorentz or gauge invariance and, respectively, 
breaking another one. This prescriptions should break the supersymmetry, 
since they do not restrict the bosonic spectrum as well. 

However, in contrast to lattice gauge models, one may not need to eliminate 
such doubling states. Since there is a conjecture that IKKT matrix model 
nonperturbatively describes both IIB and IIA string models [1], which have 
different chirality content, it may be possible that the doubling in the language 
of matrix models is related to the string duality. 
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